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Abstract. A charge motion in an electric field in a Holstein molecular chain is modeled in the absence of 
dissipation. It is shown that in a weak electric field a Holstein polaron moves uniformly experiencing small 
oscillations of its shape. These oscillations are associated with the chain's discreteness and caused by the 
presence of Peierls-Nabarro potential there. The critical value of the electric field intensity at which the 
moving polaron starts oscillating at Bloch frequency is found. It is shown that the polaron can demonstrate 
Bloch oscillations retaining its shape. It is also shown that a breathing mode of Bloch oscillations can arise. 
In all cases the polaron motion along the chain is infinite. 
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PACS. 63.20.-e Phonons in crystal lattices 
87.15.hj Transport dynamics 



87. 15. A- Theory, modeling, and computer simulation 



Presently, there are a lot of papers devoted to modeling 
the motion of a charged particle in various-type molecu- 
lar chains (see books and reviews [T]-[S] and references 
therein) . Interest in the problem is associated with possi- 
ble practical applications of such chains in nanoelectronics 
[5],[7]- Functioning of nanoelectronic devices is based on 
the motion of charges along molecular chains placed in 
an external electric field. Notwithstanding the great num- 
ber of publications on the problem of charge motion in a 
molecular chain, it would be an exaggeration to say that 
the problem is conclusively solved. 

It is well known that an electron in an ideal non- 
deformable molecular chain placed in an electric field demon- 
strates Bloch oscillations |H]-[II]- In a deformable molec- 
ular chain an electron turns into a polaron state. In the 
simplest case of a Holstein molecular chain in the absence 
of an external field the motion of an electron with steady 
velocity v appears to be impossible. The reason is that a 
steadily moving polaron state in a crystal with dispersion- 
free phonons is not the eigenstate of the Holstein Hamilto- 
nian [12]. Using a discrete Holstein molecular chain as an 
example, in |13| we showed that in a weak electric field a 
polaron as a whole demonstrates Bloch oscillations which 
in a strong electric field take the form of breathing mode of 
Bloch oscillations. In contrast to this, in [H] it was shown 
that in a continuum approximation in an electric field not 
exceeding a certain critical value, the polaron motion can 
be steady. 

The aim of this paper is to model numerically the 
charge steady motion in a discrete chain in an electric field 
and to study the conditions under which it turns into oscil- 
latory motion. In the Holstein model under consideration 
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the electron motion in a molecular chain is described by 
Hamiltonian [T4]-[l8]: 

H = i^E„(a;!:an-i + a+a„+i) + a' E„ ^nOn «" 
+ En eSana+an + E„ Pl/'^M + E„ Kql/2, (1) 

where v is the matrix element of the charge transition 
from the n-th site to the n ± 1 site, a^,an are the op- 
erators of the birth and annihilation of the charge at a 
site with number n, a' is a constant of the charge interac- 
tion with site displacements g„, e is the electron charge, 
£ is the electric field intensity, a is the lattice constant, 
P„ = —ihd/dqn, M is the site's effective mass, K is the 
elastic constant. 

To study the dynamics of the charge motion in the 
chain use is usually made of a semiclassical description 
based on the fact that the site's mass M exceeds greatly 
the electron mass. In order to pass on to semiclassical 
description we will seek the wave function of the system 
\^) in the form of its expansion in terms of coherent states: 

W{t))^hn{t)ate^v{-\Y.[PAt)^,-<t)<l3]]\^\ (2) 

n j 

where |0) is the vacuum wave function, while the quanti- 
ties /3j(t) and TTj{t) satisfy the relations: 

{W{t)\q„\^it)) = Pn{t), {^it)\Pn\^it))=7rn{t). (3) 

Dynamical equations for the quantities b„{t) and /3„(i) 
obtained from ((Ij-© have the form: 

ihbn = i^{b„-i + b„+i) + a' Pnbn + e£anb„, (4) 
M/3„ = -i<f/3„-a'|6„|2. (5) 
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Fig. 1. Motion of the polaron at a constant velocity in an elec- 
tric field of intensity E = 0.01 at uj — l,ri — 1.276 for various 
values of the parameter k: (a) k = 2, (b) k — 3. Oscillations of 
the polaron shape whose amplitude depends on the values of 
K are conspicuous (here and elsewhere r = 10~^^sec.) 



Equations (|31) are Schrodinger equations where 6„ are am- 
plitudes of the probabilities of the charge localization at 
the n-th site; ^ are classical motion equations for sites 
displacements /?„. 

Turning to coherent states using transformation ^ 
corresponds to the idea of the strong coupling polaron. 
Herewith the quantities /?„ determine the phonon "envi- 
ronment" surrounding the particle. Notice that in the ab- 
sence of an electric field, the steady motion of the polaron 
along the chain is impossible. The reason is that the group 
velocity of phonons Vg = dil/dk is equal to zero since the 
quantity f2, according to ([5]), does not depend on the wave 

vector k: fl - 



K I M . In other words, in the absence of 
dispersion in the chain, the polaron " environment" cannot 
follow the particle's motion along the chain. 

For numerical modeling of the polaron motion, we used 
dimensionless variables in which system (II]),© has the 
form: 



i— 5- = rj{b.n + l + bn-l) + KUj'^Unbn + Enbn, (6) 

at 



dt^ 



-oj'^u,-, 



(7) 



where dimensional quantities involved in (H)),® are re- 
lated to dimensionless quantities involved in (jG]),© as: 

r; = Tv/h, bj"^ = t'^K/M, q„ = /3w„, t = tI, (8) 
Kcj2 = T^{a'f/Mh, P = r'^a' /M, E = Eear/h, 

where r is an arbitrary time scale. In view of arbitrariness 
of r any set of dimensionless parameters r],LLj,K,E corre- 
spond to a continuum of possible dimensional parameters 
of the molecular chain: i', K, M, a' ,£. 
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Fig. 2. Numerical (•) comparison of the oscillation period Tpjv 
with the quantity 1/V for E = 0.02, 0.015, 0.01, 0.005, 0.003, 
0.001, 0.0007, 0.0005, 0.0003, 0.0002, 0.0001 at the following 
parameter values; n — 2, uj = 1,7] = 1.276. The right top point 
(•) on the graph corresponds to the value E = 0.0001. The 
straight line corresponds to the graph Tpjv = 1/V . 



In the continuum limit, solution of (HI),© yields the 
following relation between the polaron equilibrium veloc- 
ity V = va/r and the electric field intensity |12) : 



E = 271 



4 9 
2^ V 



1 



1 



K V^ smh^i2mj uj/kV)' 



(9) 



According to [12, the steady motion of a polaron in the 
conservative system under consideration in an electric field 
appears possible due to the fact that while the polaron 
moves along the chain it leaves behind a " tail" of excited 
oscillators. In this process the work done by the electric 
field is expended for excitation of these oscillators. 

Figure [T] shows the polaron motion in an electric field 
for various values of the parameter k, obtained as a result 
of numerical solution of discrete equations dl]),®. 

The calculations were carried out by a standard 4-th 
order Rungc-Kutta numerical method. The initial values 
of I 6n(0)| (6„ = Xn + iyn) were chosen in the form of an 
inverse hyperbolic cosine: 



|6n(0)| 

4 = i&n(o)i(-i: 




(10) 



y^ = |6„(0)|(-l)"+Vy2, 



,0 



^n — I ^n(0)P/'^ni '^"n/'^^ = 0' Thc parameter values were 
taken to be: w = 1,?7 = 1.276. The chains length was 
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fitted to suit the calculations duration with the proviso 
that the polaron should not come too close to the end 
of a chain when calculations are over. To calculate the 
numerical dependence of the polaron velocity V on the 
electric field intensity E the sequence length was taken in 
the range from N = 5000 to iV = 30000 depending on E. 
The value of tt-q (the center of an initial inverse hyperbolic 
cosine) in pUj) was fitted so that at the beginning of the 
calculations the polaron be far away from the end of the 
chain. For the parameter values chosen in this work it was 
appropriate to take uq from no = 300 to no = 700. 

To calcalculate numerically the polaron velocity V we 
used the relations: 



V = dX/di, X(i) = V n|6„(t)|^ 



(11) 



Figure [T] suggests that while moving steadily, the po- 
laron shape oscillates periodically every Tpn = l/V due 
to the presence of Peierls-Nabarro potential in a discrete 
chain. The amplitude of these oscillations turns to zero in 
the continuum approximation when the polaron size ex- 
ceeds considerably the distance between neighboring sites. 
Numerical comparison of the oscillation period Tp^{E) 
with the quantity l/V{E) for various values of the elec- 
tric field intensity E is given in Figure [2] 

Figure H^ demonstrates the dependencies of the po- 
laron velocity V on the field intensity E determined by 
relation ([9]) for various values of the parameter k. All 
the V{E) curves have an infinite derivative at zero point: 
V'(0) — oo which corresponds to infinite mobility of the 
polaron in the conservative system under consideration 
determined by Hamiltonian ([1} . Figure [3}d compares the 
average velocity of the polaron steady motion obtained 
by numerical integration of system dU,© with the theo- 
retical curve determined by (jH]) for the following param- 
eter values: k = 1,77 = 1.276, w — 1. For these values of 
the parameters n and rj, continuum approximation is ful- 
filled to a high accuracy since the polaron state radius: 
f = r/a — Arj/ K, « 5.1 ^ 1. It is seen from Figure [3)d that 
the difference between the theoretical and experimental 
values increases as V grows. This difference is likely to be 
due to the fact that in deriving the theoretical dependence 
^ we did not take into account that the polaron shape 
changes as velocity V increases. 

At small values of the electric field intensity, when the 
polaron velocity is small, its motion is not accompanied 
by Bloch oscillations Figure |4^. The reason is that the 
energy of excited oscillators in the "tail" of the moving 
polaron increases during the period of Bloch oscillations 
and exceeds greatly their energy. It can be said that while 
moving along the chain a polaron emits phonons. In a con- 
servative system, such emission plays the role of friction, 
prohibiting the polaron from accelerating in an electric 
field. This makes possible polaron's motion at a constant 
velocity in a weak electric field. 

As the electric field intensity increases, the situation 
changes: the energy of Bloch oscillations grows, while the 
length of the "tail" which has grown in the course of oscil- 
lations, decreases. When the energy of Bloch oscillations 
exceeds the "tail's" energy, they become visible FigureHJa. 
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(b) 



Fig. 3. (a) - Dependence @ of the polaron velocity V on the 
electric field intensity E for various values of the parameter 
K (k = 1,2, 3,4), r; = \.21%,bj = 1. (b) - Comparison of the 
theoretical ( — ) ^ and numerical (•) (0,(0 dependence of 
the polaron velocity V on the electric field intensity E for k = 
1,77 = 1.276,0; = 1. 



Depending on the values of the system's parameters Bloch 
oscillations can look like oscillations of a polaron as a 
whole, or may take the form of a breathing mode of Bloch 
oscillations. 

The period of Bloch oscillations in dimensionless vari- 
ables is Tbo = '^T^/E {Tbo = 2nn/T£ea, E = Sear/h). If 
the electric field intensity E = 0.05 the period of Bloch os- 
cillations is Tbo ~ 125.66. Figure |3|3 demonstrates some 
periods of Bloch oscillations at the onset of the polaron 
motion in the electric field of intensity E = 0.05. The 
graph of the function |6„(t)p on Figure [3^ demonstrates 
a small segment of the polaron movement at a constant 
speed in the electric field of intensity E = 0.0005. Fig- 
ure 0^ also represents the graph of the function X{t) pT|) 
for the time t = 34400. The length of a chain for this 
case was chosen to be iV = 32000. The period of Bloch 
oscillations for the electric field intensity E = 0.0005 is 
Tbo ~ 12566. Hence we observe a stationary motion of 
a polaron (without any oscillations) for the time corre- 
sponding to almost three periods of Bloch oscillations for 
E = 0.0005. At the end of calculations {t = 34400) we ob- 
serve the same graph of the function |6„(t)p, as presented 
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Fig. 4. Dependence of the charge motion character on the 
electric field intensity E for k. = 1,lj = l,i] = 1.276. When 
E = 0.0005 (fig. (a)) the polaron moves at a constant velocity, 
for E — 0.05 (fig.(b)), the initial polaron motion at once turns 
to Bloch oscillations and then its motion takes the form of a 
breathing mode of Bloch oscillations. 



in Figure|3^ (t = 1000). These results are in a good agree- 
ment with theoretical conclusions ^12 , Figure [3Jd. 

In order to estimate numericaUy the critical value of 
the electric field Ecr at which Bloch oscillations (BO) arise 
let us set the energy of oscillators on the interval which 
the polaron runs during the period of Bloch oscillations: 
Wosc = 'huJiT-osc (where Uosc — vTbo/ o. is the number of 
oscillating sites on the interval vTbo : Tbo is the period of 
Bloch oscillations) equal to the energy of Bloch oscillations 
Wbo = Ji^BOi ^^BO — e£a. As a result, the field's critical 
value will be Ecr = V^ttloV. The electric field is weak: E < 
Ecr, if the energy of Bloch oscillations is less than the loss 
of energy by polaron (caused by excitation of oscillations 
in the chain while the polaron moves along it) during the 



period of the polaron's BO. Accordingly, the field is strong 

if ^ > Ecr- 

The upper branch of the curve for the V{E) depen- 
dence in Figure [5^ is unstable and is not reproduced in 
experiments. When E > Emax, where Emax is the maxi- 
mum electric field at which the steady motion of the po- 
laron is possible, (in numerical experiments this happens 
well before E = Emax is achieved) Bloch oscillations be- 
come dominating. In this case the energy which the elec- 
tron gets from the electric field goes into Bloch oscillations 
which, in turn, transfer it to the oscillators of the chain. 

The velocity at which the mass center of a Bloch oscil- 
lator moves in the electric field is also limited and reaches 



its maximum at E < E„ 



In the limit i? — >■ oo this ve- 



locity vanishes, since the BO amplitude turns into zero 
and the polaron localizes at one site. 

The results obtained testify to the complicated char- 
acter of the charge motion in a molecular chain. 

The statement that a Holstein polaron cannot move in 
a molecular chain with dispersion-free oscillations of sites 
was first made by Davydov and Enol'skii in [T^. In their 
subsequent works this statement was generalized to the 
case of polar crystals with optical phonons [20], [21]. Some 
erroneous works devoted to this problem were analyzed in 

In conclusion it may be said that in this work we have 
carried out a direct numerical modeling of the Holstein 
polaron motion in a dispersion-free chain. The modeling 
demonstrated that in an electric field E < Ecr the po- 
laron moves along the chain at a constant velocity. In the 
range of parameter values at which the polaron size be- 
comes comparable with the lattice constant, its motion is 
accompanied by oscillations of its shape (height and width 
of distribution |&nP)- In the course of oscillations the po- 
laron velocity V remains invariant. In the conservative 
system under consideration this invariance is caused by 
the fact that while traveling along the chain, the polaron 
excites oscillations of sites. 

When E > Ecr , the steady motion of the polaron be- 
comes impossible: the energy obtained by the polaron from 
the field per unit time becomes greater than that lost to 
the excitation of oscillations in the chain. The steady mo- 
tion in this case "falls" and turns to Bloch oscillations. 
As distinct from Bloch oscillations in a non-deformable 
chain when the oscillation center is left in place (finite 
motion), in a deformable chain the polaron's mass cen- 
ter demonstrates translational movement accompanied by 
Bloch oscillations (infinite motion). In a strong field the 
polaron motion becomes similar to the breathing mode of 
Bloch oscillations. 

Though in this work we have considered the case of 
a one-dimensional Holstein chain, the results obtained are 
rather general and can be applied to any media with dispersion- 
free phonons in which formation of a localized polaron 
state is possible. In particular, these results are valid for 
a strong coupling polaron in ionic crystals. 
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